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psj I Abstract 

(— I I A relic lepton asymmetry orders of magnitude bigger than the baryon one may hide in the rehc 

Qh' neutrino background. No direct theoretical or experimental limitations on its magnitude and sign 

are known. Indirect cosmological constraints exist ranging from \L\ < 0.01 to L < 10. We discuss a 
D I BBN model with late electron-sterile neutrino oscillations which is a fine leptometer - it is capable 

' of feeling extremely small relic lepton asymmetry - |L| > 10^^. This sensitivity is achieved through 

the influence of such small L on the neutrino oscillations, suppressing or enhancing them, and thus 
changing the primordially produced ^He. The influence of L on nucleons freezing in pre-BBN 
I epoch is numerically analyzed in the full range of the oscillation parameters of the model and for 

!>■ ■ L> IQ-io. The case of oscillations generated asymmetry by late electron-sterile oscillations and 

its effect on primordial '^He is also briefly discussed. 

1 Introduction 

Lepton asymmetry of the Universe is usually defined as L = {Ni — Ni)/N^, where Ni is the number 
^ ' density of leptons, and Nj of antileptons, while is the number density of photons. In case of 
c3 equilibrium ^ = /x/T, the degeneracy parameter, is used instead as a qualitative measure of the lepton 
asymmetry. 

It is traditionally assumed that the lepton asymmetry is of the order of the baryon one, which is 
measured precisely by different independent means (BBN, CMB data) to be /3 = {Nb—Ni)/N^ ~ 6.10~^°. 
However, actually big L may reside in the neutrino sector and hence, L may be many orders of 
magnitude larger than /3. Therefore, L is defined mainly by the sum of the asymmetries in the different 
neutrino sectors L ~ J2Lu,- Direct measurements of the lepton asymmetry magnitude and sign have 
not been done yet. There are available just indirect indications and constraints. 

Studying L of the Universe is intriguing and important because of many reasons. Just a short list 
of some of them reads: 

• A significant cosmological effect of L might be expected, having in mind that neutrinos are 
abundant and large L may be contained in the neutrino sector. 



^Universal charge neutrality implies that the lepton asymmetry in the electron sector is Le ~ /3. 
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• Precise determination of neutrino properties is of cosmological importance. As far as the uncer- 
tainty of neutrino characteristics leads to large systematic errors in the estimation of cosmological 
parameters obtained from CMB data. Therefore, knowledge about L and its nature would help 
to determine them more precisely. 

• Studying L and its cosmological influence provides, on the other hand, an opportunity to use 
cosmology as a probe of neutrino properties. 

• There exist different mechanisms of generation of L, and among them is the natural possibil- 
ity of amplifying the neutrino asymmetry by neutrino active-sterile resonant oscillations in CP- 
asymmetric plasma of the early Universe [II |2]. 

• Planck will soon be able to better constrain L, namely Planck will measure the radiation content 
at CMB decoupling with a precision SNf,ff = ±0.26, i.e. it will be able to measure with this 
precision the eventual L contribution to the radiation content. Besides, in case ^ > 2 Plank will 
be able to detect small neutrino mass of the order of ~ 0.07 eV because L enhances the effect of 
the mass. 

• Today it is reasonable to study L, also as a possible solution of the recently found cosmological 
preference or /and indication for additional relativistic density [3]. Namely, recent measurements 
of primordially produced ^He Yp |3] and CMB measurements of Yp [5J point to an effective number 
of the relativistic degrees of freedom at the BBN epoch Nbbn = 3.8lo 7, at the CMB formation 
epoch NcMB = 4.341^:88 at 68%, and on the basis of LSS survey data Nlss = 4.8111 at 95%, i.e. 
higher values than previously estimated ones and higher than the theoretically predicted standard 
value for 3 neutrino species N^jf = 3.046. 

• Moreover, recent analysis of the combined neutrino oscillation data, including LSND and Mini- 
Boone requires 1 or 2 additional low mass sterile neutrinos z/^, participating into oscillations with 
flavor neutrinos with higher mass differences values, than the ones required by solar and atmo- 
spheric neutrino oscillations experiments. These light z/^ brought into equilibrium by oscillations 
with active neutrinos may be a successful explanation of the excess relativistic density. Besides, 
it is known that such active-sterile oscillations in the resonant case may generate L in the CP- 
asymmetric plasma of the Universe. 

• Knowledge about L magnitude and sign is relevant for cosmological models with leptonic domains, 
for inhomogeneous BBN models, baryogenesis through leptogenesis issues, etc. 

• Determining L at BBN epoch would allow testing the assumption that sphalerons equilibrate 
lepton and baryon asymmetries. 

• Particularly concerning neutrino oscillations: Due to the fact that L is capable of suppressing and 
inhibiting or enhancing neutrino oscillations, determining L will enlighten our knowledge about 
their cosmological role. 



In the next section lepton asymmetry effects on processes in the Universe and some of the based 
on them cosmological constraints on L are shortly discussed. In section 3 the interplay between small 
lepton asymmetry and electron-sterile neutrino oscillations u^, Ug, effective after neutrino decoupling, 
is described in terms of general kinetic equations for oscillating neutrinos and numerically analyzed. In 
section 4 the results of our numerical analysis of the influence of relic L and oscillations generated L on 
primordial production of ^He in case of late o Us oscillations are presented. 
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2 Knowledge about Lepton Asymmetry 



Cosmic neutrino Background has not been detected yet, hence L is measured/constrained only indi- 
rectly through its effect on other processes, which have left observable traces in the Universe, like the 
abundances of light elements produced in BBN, Cosmic Microwave Background, LSS, etc. (see for 
example the reviews [6] - [9]). Recently it was found possible to obtain an information about L from 
QCD transition, on the basis of L effect at QCD epoch [TB]. Below we review the main cosmological 
effects of L and some of the recent cosmological constraints on L magnitude and sign. 

2.1 Lepton Asymmetry Effects 

A. Energy density increase. 

A well-known cosmological effect of L is the increase of the radiation energy density, which is usually 
expressed in terms of the increase of the effective number of the relativistic degrees of freedom 

where and are the photon and neutrino energy densities, correspondingly. In equilibrium L may 
be expressed as usual through the chemical potential p or degeneracy parameter ^ = p/T: 

i 

The increase of N^ff due to L is 

The increase of the radiation density due to L speeds up the Universe expansion H = (S/SvrGp)^/^, 
delays matter/radiation equality epoch, changes the decoupling temperature of neutrino, which on their 
turn influence BBN, CMB and the evolution of the density perturbations, i.e. formation of structures 
in the Universe. 

Particularly well studied is L effect on BBN. The increase of the cooling rate of the Universe due to 
L, leads to earlier freezing of the reactions governing neutron-to-proton ratio n/p , i.e. leads to higher 
freezing ratio {n/p)f, which reflects in higher D and '^He abundances. 

B. Direct kinetic effect. 

Besides its dynamical effect, lepton asymmetry with a magnitude \L\ > 0.01 in the Uf, sector exerts 
also a direct kinetic effect on the n-p kinetics and on BBN, because the Uf, participates in the reactions 
interconverting neutrons and protons. In this case the effect on BBN and the outcome of the light 
elements is L sign dependent. 

As is obvious, L > in the pre-BBN epoch would result into reduction of {n/p)f and thus leads to 
light element underproduction, while L < would lead to their overproduction. Degenerate BBN has 
been thoroughly studied (see for example pioneer papers ref. |13j). 

An empirical formula, which provides a fairly good fit (see for example ref. [Hj), presents the 
dependence of the produced primordially *^He, Yp, on the discussed dynamical and kinetic effect of L: 

Yp ~ (0.2482 ± 0.0006) + 0.0016r/io + OMSAN^fj - 0.3^^^ 

C. Indirect kinetic effect due to asymmetry-oscillations interplay 

Small L, \L\ << 0.01, that has negligible A and B effects, in case of late z/g ^ oscillations may 
considerably influence oscillating Ue, namely change its evolution, number density, energy distribution, 
oscillation pattern and thus through Ue influence BBN kinetics [21 [17]. The effect of small relic L and 
nonresonant i^s oscillations effective after neutrino decoupling on BBN has been first studied in 
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ref. [IT]. It was found that L < 10^'' is destroyed by oscillations, while L > 10^'^ may enhance or 
suppress oscillations and through them influence primordially produced elements (see also ref. 

On the other hand active-sterile oscillations may induce neutrino- antineutrino asymmetry growth 
during the resonant transfer of neutrinos [H |21 [iHl [19]. The case of weak mixing and relatively big 
mass differences > 10~^ eV^ was discussed first in ref. [1], while the case of asymmetry generated 
at relatively big mixing and small mass differences (5m^ < 10~^ eV^ was first found in ref. |2]. This 
dynamically produced asymmetry exerts back effect on oscillating neutrino and changes its oscillation 
pattern. □ There exists another possibility: when L growth is not high enough to have a direct 
L kinetic effect on the synthesis of light elements, it can effect indirectly BBN through its effect on 
oscillating neutrinos. In the case of resonant ^ Vg oscillation effective after neutrino decoupling 
neutrino-antineutrino asymmetry is amplified by not more than 5 orders of magnitude from an initial 
value of the order of the baryon one |T9j. Oscillations generated asymmetry suppresses oscillations at 
small mixing angles, leading to noticeable decrease of ^ife production at these mixing angles. The effect 
of small L generated by oscillations on ^ife abundance and on cosmological constraints on oscillations 
was analyzed in ref. [2l |19l |20] . 

In conclusion, very small asymmetries 10~^ < L <^ 0.01 , either relic or produced in active-sterile 
oscillations, may considerably influence oscillating electron neutrino and through it Yp and BBN. In 
the next sections we study numerically this special case of influence of small L on oscillating neutrino 
and on BBN and show that BBN produced '^He feels extremely small L and represents now the finest 
known " leptometer" . For comparison, in the next subsection we present first the available cosmological 
constraints on L coming from BBN with flavor oscillations present. 

2.2 Lepton Asymmetry Constraints 

At present BBN provides the most stringent constraints on L. There exist numerous papers on the 
subject. I list here only the results of more recent work. For more information and reference of earlier 
papers see the review paper [B]. 

In case of equilibration of the neutrino oscillation degeneracies due to flavor oscillations before BBN 
the limit on L in the muon and tau neutrino sector strengthens. Then BBN constraint reads [23] 
\S,v\ < 0.1. A more recent study for that case gives the constraint : —0.04 < ^^j^ < 0.07 [24J. For such 
a small L the expantion rate remains practically standard because SN^ff 10^^. 

However, it was realized later that the equilibration of the chemical potentials before BBN depends 
on the value of the yet unknown mixing ^13. [25j. Besides, since relic L is capable to suppress or 
enhance oscillations, depending on its value and the values of oscillation parameters, L itself may play 
the role of an inhibitor or a catalyzer of the equilibration. Hence, different possibilities for the chemical 
potential in different neutrino flavors still may have place. 

The analysis on the basis of BBN and D and '^He abundance and CMB/LSS constraints on baryon- 
to photon value, provided restrictive constraints on the neutrino degeneracy [Hj. Namely the following 
constraints were derived for iVe// = S.Stg;^ and different possibilities for the chemical potentials: in case 
^lye 7^ = iv < 2.3 corresponding to L < 5; in case = 7^ iv^ < ^ L < 7.6, while in case 

= iu^, = C,u-r 0.01 < < 0.1 and L < 0.07. In the last case practically the rate of expansion does not 
change, and the small dynamical effect of L corresponding to AN^ff ~ 0.03 is undetectable by BBN 
and CMB [14j. 

CMB and LSS provide much looser bounds: L modifles the power spectra of radiation and matter 
[9], [15]. However, today's sensitivity of CMB and LSS data does not allow to probe different flavors. 

^ The possibility of L to suppress neutrino oscillations was discussed in the case of oscillations with larger mass 
differences first in ref. |21) . 

•^The case of fast oscillations generating high enough L that exerts direct kinetic effect B on BBN was discussed in 
many publications (see for example ref. [22)). 
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They feel only the change in the total density, i.e. SN^ff. CMB data (including WMAP 5 years data 
results) combined with LSS data puts the bound: < 0.7 and L < 0.6 at 2a level. 

The WMAP5 data combined with the data on BBN produced helium-4 provides more stringent 
bounds, namely: —0.04 < < 0.02 in case of equilibration, while otherwise —0.03 < < 0.13, 
I^i'mtI < 1-67 [lOj. See however the recent work [26] based on the new data on Yp and WMAP7, which 
allows considerable relaxation of the constraints, namely —0.14 < < 0.12. 

In conclusion, depending on the different combinations of observational data sets used and the 
assumed uncertainties, cosmology provides an upper bound for L in the range |-^^z/p,rl < 10~^ — 10 and 
\L\ < 0.01 — 0.2. These values are many orders of magnitude larger than the baryon asymmetry value. 

In the next sections we analyze the effect of small L on BBN with late neutrino oscillations and 
present a possibility for much sensitive leptometer capable to constrain L values closer to the baryon 
one. Namely, we describe a model of BBN with late electron-sterile neutrino oscillations, which can 
'measure' L as small as L = 10~®. 

3 Small Lepton Asymmetry and Late Neutrino Oscillations 

Here we present the results on the interplay between small asymmetries L << 0.01 and late Ue ^ Vs 
oscillations in the early Universe and their effect on BBN. This study continues and broadens the 
analysis of refs. [T7] and [TU] . 

The effect of large L either previously existing or produced by oscillations, which effect directly the 
kinetic of BBN has been studied in numerous papers. The effect of small relic asymmetry was not so 
thoroughly studied. Its effect on primordial ^He abundance was first analyzed for hundreds of — 9 
combinations in refs. [17] and [21 119] . 

We have analyzed electron-sterile neutrino oscillations 

vi = cos(6')z/e -|- sin(^)z/s 

z/2 = — sm{9)h'e + cos{9)h's, 

where 9 is the mixing angle, ui and 1^2 are the Majorana particles with masses mi and 7712- The sterile 
neutrino Ug is not having the usual weak interactions, and is assumed to have decoupled much earlier than 
the flavor neutrinos, hence its density is much lower than the density of electron neutrino n^^ << Uy^. 
The case of neutrino oscillations effective after active neutrino decoupling 5rn? sin^ 29 < 10^'' eV^ is 
studied in detail. 

For that specific case of small L both the dynamical effect A, discussed in the previous section, and 
the direct kinetic effect B on nucleon kinetics are negligible. We will discuss the asymmetry effect C on 
oscillating neutrinos and through them on BBN. 

In the case of late oscillating active-sterile neutrinos with relic or generated in oscillations L the 
energy distribution of neutrinos may be strongly distorted from the equilibrium Fermi-Dirac form [19] . 
Hence, a precise account for the energy spectrum distortion of the degenerate oscillating neutrinos is 
necessary to reveal the effect of small lepton asymmetry. Particularly the capability of small relic L 
to enhance oscillations has essentially spectral character [17] and requires a precise kinetic approach, 
provided in the numerical analysis, described below. 

We have studied two different cases of L, namely initially present at the neutrino decoupling epoch, 
called further on relic L, and dynamically generated L by oscillations. 

3.1 Evolution of Oscillating Neutrino in Presence of Lepton Asymmetry 

We have used a self consistent numerical analysis of the kinetics of the oscillating neutrinos, the nucleons 
freeze-out and the asymmetry evolution for the analysis of lepton asymmetry role in BBN. Kinetic 
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equations for neutrino density matrix and neutron number densities in momentum space are used to 
describe the evolution of the system of oscillating neutrinos in the high temperature Universe, following 
the approach of ref . . 



dp{t)/dt = Hp, {dp{t)/dp,) + 

+1 [Ho, p{t)] + tV2GF {C - Q/M^) [a, pit)] + O (gI) , 
dpit)/dt = Hp, {d pit) /dp,) + 

+1 [Ho, pit)] + tV2GF (-£ - Q/M^) [a, pit)] + O , 

drin/dt = Hpn idnn/dpn) + 

+ / dfi(e",p, v)\Aie~p h'n)\^ K-?^p(l - Pll) - nnPiii^ - n^-)] 



aij = U*JJje, Vi = UiiUiil = e,s). Ho is the free neutrino Hamiltonian. Q arises as an W/Z propagator 



efrect,g r^E,T. 2L,^ + L,^ + L,^, L^,, ~ (iV^,, - iV^,,)/iV^ L,^ ~ / d^pipLL - Pll)/N. 



The first two equations describe the evolution of neutrino and antineutrino ensembles. They provide 
a simultaneous account of the different competing processes: expansion (first term), neutrino oscillations 
(second term), neutrino forward scattering and weak interaction processes. The number densities of 
nucleons and electron neutrino were assumed the equilibrium ones. The sterile state was assumed empty 
at the time of decoupling of the electron neutrino. □ 

Due to the non-zero L term the equations are coupled integro-differential and the numerical task is 
much complicated than in the case of zero L. Besides, L term leads to different evolution of neutrino 
and antineutrino due to the different sign with which it enters their equations. (The case of L P 
corresponds to negligible C term in the potential when the evolution of the neutrino and antineutrino 
density matrices is identical. 

We studied numerically the evolution of neutrino ensembles, evolution of L, and also the evolution 
of nucleons and L role during pre-BBN epoch for a broad range of oscillation parameters and 10^^*^ < 
L < 0.01 for the case of a relic L. For the case of oscillations generated L its initial value was 
assumed L 13. The numerical analysis was provided for the temperature range [0.3 MeV, 2 MeV] 
and the full set of oscillation parameters of the electron-sterile oscillations model, and with higher 
accuracy than in previous studies. We have calculated precisely neutron to nucleons freezing ratio 
Xl = nl/iun + UpY = fi5rn?,svo?26,L) which is essentially influenced by oscillations and L. The 
primordially produced ^ife was estimated from it. 

Active-sterile oscillations proceeding after neutrino decoupling produce Pg at the expense of active 
neutrino and thus 6Neff does not change. However in that case late oscillations for a wide range of 
values of oscillations parameters and L, strongly distort neutrino energy spectrum. Therefore, a precise 
description of neutrino momenta distribution is necessary. In this analysis we have used between 1000 
and 5000 bins to describe neutrino spectrum distribution in the non-resonant neutrino oscillations case, 
and up to 10 000 in the resonant case. Depending on the oscillation parameters and L values, the 
following interplay between L and oscillations can be observed: relatively large L suppress oscillations, 
smaller L lead to their resonant enhancement. On the other hand, resonant oscillations also are capable 

* the case of non-zero population of the sterile neutrino state was considered in detail for L = 10"^" case in refs. [271 




dQie^,p, iy)\Aie^n pi>)\^ K+n„(l - pll) - rippiLi^ - ne+)] . 
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to amplify L. In the following section we present the results of a detail numerical study of this interplay 
on BBN produced ^ife. 

4 Lepton Asymetry, Neutrino Oscillations and BBN. The Re- 
sults 

To study the lepton asymmetry effect on BBN we have provided a detail numerical analysis of the 
influence of L on Yp, because primordially produced ^ife is highly sensitive to the nucleons kinetics 
during the pre-BBN epoch and besides, it is the most precisely measured element among light elements 
synthesized during BBN. A recent measurement of Yp was provided on the basis of 93 spectra of 86 low 
redshift HII regions 

4.1 Oscillations Generated Lepton Asymmetry and BBN 

In the analyzed oscillations case the evolution of L is dominated by neutrino oscillations and typically 
L has rapid oscillatory behavior: it oscillates and changes sign. We have determined numerically the 
region of parameter space for which noticeable generation of LA is possible. A good approximation 
is |(5m^| sin^ 26* < 10~^'^ eV^. The maximal possible growth of L is by 4 -5 orders of magnitude. 
The instability region and the magnitude of L are close to the bounds existing in literature for other 
oscillation models [30] . 

Precise description of the distribution of the neutrino momenta was found extremely important for 
the correct determination of L evolution in the resonant oscillations case. In some cases increasing 
the resolution of momentum space leads to changes of the oscillatory character of L and diminishes L 
amplitude, This observation is in accordance with the studies of other authors in other parameter 
regions corresponding to smaller 6 and bigger [21]. Further analysis is required to decide if the 
oscillatory behavior and strong asymmetry growth is induced by numerical error. This observation 
revives the puzzle: Is the asymmetry growth due to lack of numerical accuracy?! 

L changes energy spectrum distribution and the number densities of electron neutrinos from standard 
BBN case. This influences the kinetics of nucleons during BBN and changes the production of light 
elements. 

We have precisely followed the evolution of nucleons in the presence of electron-sterile neutrino 
oscillations in the pre-BBN period for different sets of oscillation parameters and different values of L. 
The production of '^He was numerically calculated and compared to the BBN value without asymmetry 
growth account. 

Fig. [1] illustrates the typical behavior of the frozen neutron number density relative to nucleons 
when increasing the mixing in case of asymmetry growth (red curves) and in case without asymmetry 
growth account for two different mass differences. The asymmetry growth takes place at smaller mixing 
angles when increasing 5m^. Then due to the asymmetry growth the production of (correspondingly 
Yp) decreases at small mixing. The effect of the asymmetry growth on helium production is always 
towards decreasing of the caused by oscillations overproduction of ^i7e, leading to a relaxation of BBN 
constraints at small mixings. 

In case of the small lepton asymmetry values, discussed here, BBN constraints on neutrino oscilla- 
tions may noteably change [T9l [28| [29]. They become less stringent at small mixing angles (where the 
growth of asymmetry takes place). 

Thus, the analysis has proven that BBN is sensitive to the oscillations generated asymmetry, which 
usually grow not more than 5 orders of magnitude and are small \L\ < 10^^. 

^As a rule in these cases the evolution of the neutrino ensembles is strongly distorted from the expected behavior. 
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Figure 1: The dependence of the frozen neutron number density relative to nucleons on 
the mixing in case of the account of asymmetry growth (red curves) and in case without 
asymmetry growth account for two different mass differences 5im? = 10~^ eV^ and = 
10-^ eV^. 



4.2 Initial asymmetry, oscillations and BBN 

Numerical analysis of Yp{6m'^, 6, L) dependence has been provided for the entire range of mixing pa- 
rameters of the model and relic L > 10~^°. Small L, 10~^ < L << 0.01, that do not effect directly 
BBN kinetics, influence indirectly BBN via oscillations in agreement with previous analysis [17] . 

The calculated ^ife production dependence on oscillation parameters and on L shows that, in 
case of neutrino oscillations: i) BBN can feel extremely small L: down to 10"*^. ii) Large enough L 
change primordial production of '^He by enhancing or suppressing oscillations. Depending on oscillation 
values L > 10"'' may enhance oscillations, while L > 0.1(5m^/eV^)^/^ may suppress oscillations, and 
asymmetries as big as L > ((5m^/eV^)^/^ inhibit oscillations. L enhancing oscillations leads to a higher 
production of Yp. L suppressing oscillations decreases Yp overproduction by oscillations. L bigger 
than 10~^ leads to a total suppression of oscillations, i.e. to the standard BBN yield of Yp, without 
oscillations. 

In Fig. 2 the dependences of on relic L for different mixings (to the left) and different mass 
differences are presented. For L smaller than ~ 10~^ keeps unchanged from the case without 
L. The higher the mixing - the higher is the overproduction of He-4 due to oscillations. Increasing 
further L for fixed oscillation parameters leads first to an increase of helium production, corresponding 
to the region of parameters space where L enhances oscillations, and then to a decrease of helium 
production, corresponding to big L suppressing oscillations, and hence to less Yp overproduction caused 
by oscillations. At some critical L value defined by the concrete set of oscillation parameters Lc^Sm"^, 6) 
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Figure 2: Frozen neutron number density relative to nucleons as a function of the relic initial 
lepton asymmetry for 5m? = 10~^ eV^. The solid curve corresponds to maximal mixing, the 
dashed curve to sin^ 29 = lo^^ *^^ and the dotted curve to sin^29 = 10^°'^. 



the produced helium reaches its standard BBN value - i.e. L has stopped the oscillations. As is 
illustrated in the figure, the width of the enhancement region and the height of the overproduction peak 
is sensitive to the mixing. Bigger values (up to about an order of magnitude)for Lc are necessary to 
inhibit oscillations when decreasing the mixing. 

Fig. 3 illustrates the dependence of X„ on the initial asymmetry value for a fixed mixing, namely 
sin^ 29 = 10~° °^ and different mass differences Sm"^ = 10^^ eV^ and = 10~^ eV^. The enhancement 
peak due to L is more clearly expressed for higher mass differences, and is bigger for bigger mass 
differences. L bigger than ~ 10"'^ leads to a total suppression of oscillations effect on BBN for late 
oscillations studied here and hence, eliminates the BBN bounds on oscillation parameters. In that case 
instead the following approximate bound holds: 5m? /eV"^ < Li^^"^ . 

Depending on its value, relic L may also change BBN bounds: It relaxes them at large mixings and 
strengthens them at small mixings, as illustrated in ref. [17] 

The next Fig.4 presents the dependence of on the mass differences at a fixed non maximal mix- 
ing angle and for two different initial L. As illustrated, higher initial L leads to an increase of helium 
production at bigger mass differences, and reduces helium production at smaller mass differences. Cor- 
respondingly, increasing L at a fixed mixing leads to relaxation of the bounds at small mass differences 
and strengthens them fat big mass differences. 

At maximal mixing, however, bigger L leads to a suppression of the production of helium for all 
mass differences, and L = 10~^ is enough to eliminate oscillations effect, i.e. to eliminate also the 
constraints on oscillation parameters in the discussed BBN model (see Fig. 5). 

Finally in Fig. 6 we present the dependence of the helium production on the mixing angle at different 
initial L. Bigger L leads to decreasing the production of ^He with increasing the mixing. I.e. for fixed 
mass differences, L relaxes the BBN constraints at large mixings. Analysis at bigger mass differences 
5m^ > 10~^ eV^ has shown that at fixed 5m'^ L strengthens the constraints at small mixing. The results 
are in agreement with the conclusions of ref. [17], where 

the change of BBN bounds on neutrino oscillations in the presence of relic L was studied, namely: 
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Figure 3: Frozen neutron number density relative to nucleons dependence on the initial 
asymmetry for sin^ 26 = io~°-°^ and two different mass differences Sm"^ = 10~^ eV^ (lower 
curve) and = 10^^ eV^ (upper curve). 



L relaxes BBN bounds at large mixings and strengthens them at small mixings. 

In conclusion, both small asymmetry generated by neutrino oscillations and small relic asymmetry 
influences the model of BBN with oscillations, because the produced primordially '^He in this model 
feels extremely small L, namely 10^^ > L << 0.01. Hence, BBN with oscillations presents a precise 
leptometer. 

5 Summary 

The lepton asymmetry of the Universe may be much bigger than the baryon asymmetry, and hidden 
in the neutrino sector. Since relic neutrino background is not yet detected, the lepton asymmetry in 
the neutrino sector may be measured/constrained just indirectly, namely by its influence on Universe 
expansion. Big Bang Nucleosynthesis, Cosmic Microwave Background, LSS, etc. 

We discuss the case of small lepton asymmetry influence on the neutrino involved processes in 
pre-BBN epoch, and particularly on neutrino oscillations and BBN. 

We have provided a detail numerical analysis of the interplay between small lepton asymmetry 
L << 0.01, either relic (initially present) or dynamical (generated by MSW active-sterile neutrino 
oscillations) and oscillating neutrino for the case of active-sterile neutrino oscillations occurring after 
electron neutrino decoupling. The evolution of asymmetry growth in case of small mass differences 
and relatively big mixing angles was studied in more detail. Higher resolution for the description 
of the neutrino momenta distribution is required for the investigation of the asymmetry behavior in 
this oscillation parameter region. The instability region in the oscillation parameter space, where 
considerable growth of L takes place, was determined numerically. In the case of relic lepton asymmetry 
we have determined the parameter range for which L is able to enhance, suppress or inhibit neutrino 
oscillations. 

Cosmological influence of such small lepton asymmetries, which do not have direct effect on nucleons 
kinetics during BBN and are invisible by CMB, is discussed and shown to be considerable. Lepton 
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Figure 4: The dependence of the frozen neutron number density relative to nucleons on the 
mass differences at sin^ 26 = 10~°'^ and for two different initial lepton asymmetries L = 10~^° 
(the dashed curve) and L = 10"^ (solid curve). 



asymmetries as small as 10~^ may be felt by BBN in case of neutrino oscillations. The effect of 
the dynamically generated and initially present L on BBN with oscillations was studied. Relic L 
present during BBN, depending on its value, may increase, decrease overproduction of Yp or reduce it 
to the standard BBN value. Correspondingly, it can strengthen, relax or wave out BBN constraints 
on oscillations. It relaxes BBN bounds at large mixing and strengthens them at small mixings. Large 
enough L alleviates BBN constraints on oscillation parameters. In that case, instead, L constraint on 
oscillation parameters are derived. 

Oscillations generated asymmetry at small mixing angles decreases the production of Yp and relaxes 
BBN constraints at these angles. 

The discussed model of BBN with late neutrino oscillations is sensitive to extremely small L. 
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